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Abstract 

We study the effect of dipolar interactions on a magnetic striped monolayer with 
a microscopic unit cell of square symmetry, and of size (N x x N y ) spins. Even if the 
aspect ratio r = N x /N y is very large, an in-plane shape anisotropy is always negligible, 
except if N y is fairly small (N y < 40). In-plane domains are not possible, except for 
values of the dipolar coupling larger than the domain wall energy. 

1 Introduction 

An important and well known feature of magnetic dipole-dipole interactions is that, in spite 
of their weakness with respect to exchange coupling, they play an important role in magnetic 
systems [|IJ. In fact, the long range character of the magnetostatic interactions is relevant 
for determining both the ground state of the system and the excitation spectrum @; 
perhaps, the most striking consequence of dipolar interaction is the breaking of a bulk 
sample in several magnetic domains. 

The peculiarity of dipolar interactions in three dimensions is elucidated by the fact that a 
shape anisotropy is always present, independently of the size of the sample: as a consequence 
of this, it has no meaning to speak about an "infinite sample" without specifing its limiting 
shape. For example, for the slab geometry (see Fig. lb), the demagnetizing factors are in 
the x, y directions and 1 in the z one, perpendicular to the slab. Thus, the in-plane collinear 
state has a lower energy than the perpendicular collinear one, whose surface energy density 
turns out to be proportional to the thickness of the slab, and therefore the shape anisotropy 
per spin of a slab is a constant and doesn't depend on its thickness N z \ 

In genuine two dimensional magnetic systems, like magnetic films of atomic thickness 
(Fe, Co, Ni, Gd) grown on nonmagnetic substrates (Au, Ag, Cu, . . .), the situation is quite 
different: the shape anisotropy due to magnetic dipole-dipole interactions favours an in- 
plane magnetization without establishing any preferential direction within the film H|, 
and the magnetization direction of a collinear configuration is determined by the competition 
with possible easy-axis anisotropies generated by the breaking of the translational invari- 
ance in the growth direction and which may favour a perpendicular state. For in-plane 



magnetization, a monodomain configuration is energetically favoured; for perpendicular 
magnetization, the system prefers to break into domains || [7j], however weak the dipolar 
interaction may be. 

Recently, two-dimensional mesoscopic structures || have attracted attention for their 
potential technological applications. They are characterized (see Fig. la) by a "small" lateral 
dimension in one direction (N y <C N x : magnetic wires) or in both directions (magnetic dots). 
Such structures will be called striped monolayers in the following, and N x , N y will be allowed 
to assume any value. 

The aim of the present paper is to elucidate the role played by magnetic dipole-dipole 
interactions in such systems, provided that the magnetization lies in the plane of the stripe 
and the microscopic unit cell has square symmetry, and to ascertain whether the results 
valid for a slab apply analogously, namely: (i) the presence of a strong shape anisotropy 
favouring the x direction if N x ^> N y , and (ii) the existence of domains, if some in-plane 
magnetocrystalline anisotropy forces the magnetization along the y direction. As we will 
show, neither of the previous points holds. In contrast, we find that (i) even for a stripe 
with an infinite aspect ratio (N x /N y — > oo), the shape anisotropy per spin vanishes upon 
increasing N y , and (ii) in-plane domain structures can appear only for very large values of 
the dipolar interaction Q. 



2 Shape anisotropy of a striped monolayer 

Let us start by writing down the dipolar energy for an in-plane collinear ground state: 



^dip 2 2^ 3 

n,m nm 



^'nm 



(1) 



where the distance r nm = \n — rh\ (with n ^ m) is measured in units of the square lattice 
constant. 

The shape anisotropy is determined solely by the anisotropic part of Ea p (the second 
term in square brackets), which becomes, if S = (cos 6, sin 8) and up to a constant: 



E dip = -30 cos 2 5 V " mJ = -30 cos 2 9-S, (2) 

n,m r nm 

where the "cross" term of the form J2n,m( x nmynm/r nm ) vanishes if the system has at least 
one symmetry axis (a rectangular monolayer has two). 

The previous expression vanishes for an infinite two dimensional monolayer, or for any 
square-shaped system (N x = N y ), provided that the symmetry of the microscopic unit cell 
is square. For a rectangle with N x > N y , a couple of points (n, m) will generally correspond 
to a value \x nm \ > \y n m\, so that the summation S will have a strictly positive value, E dip 
will be minimal for 8 = and the magnetization will be oriented along the x axis. 

The relevant quantity is the dipolar energy per spin: e^p = Edi P /N x N y , and we are 
interested in analyzing the dependence of on N x and N y . For N y = 1, in the limit of an 
infinite chain of spins, it is immediately found that 

t dip {N x = oo, N y = 1) = -30 cos 2 . i £ = -3fi cos 2 6 ■ C(3) , (3) 
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where £(3) = 1.2021 is the Riemann's zeta function ||. 

The competition between the dipolar contribution (3£(3)fl) and possible in-plane ani- 
sotropies favouring the y axis, determines the actual direction of the magnetization. In the 
following, we will evaluate £<&p in a continuum approximation, for any value of N x , N y . 

The exact definition of S is: 

S = - 2 Y. E 77J ' / with (4) 

1 l,l'=l m,m'=l [{I ~ I f + vn- m'Y\ 

It is useful to introduce the new integer variables: x = (I — I'), with — (N x — 1) < x < 
(N x — 1) and y = (m — m'), with — {N y — 1) < ?/ < (iV y — 1). 5 now writes as a sum on 
solely two indices: 

1 JV„-1 N y -1 2 _ 2 

5 = 2 E E (iv,-N)(iv,-H) J (5) 
z x=-(N w -i)y=-(N y -i) y-t ^ y ) 

At this point we exploit the fact that S = if N x = N y . Therefore, we will write 
N x = N y + (N x - N y ): 

s = ~ E (*»-M)(tfv-lvD- 



1 |x|<(JV«-l) 2 _ „2 

E E (N x -\x\)(N y -\y\y ■ ' 



z ix|>(jv v -i) iwi<(jv„-i) ^ ^ </ ; 

+ l(N x -N y ) (N y -\y\)- J ' 1 ~" 1 



Z \x\,\y\<(N y -l) \ x +V ) 

The first sum vanishes, as seen by interchanging the two dumb indices (x, y); in the third 
sum, for the same reason we can get rid of the quantity proportional to (N x — N y )N y (x 2 — 
y 2 )/(x 2 + y 2 ) 5 / 2 , whilst in the other one the term y = does not contribute. Finally, by 
using the parity of the addenda, we obtain: 

Nx -l N y -1 2 _ ..2 

^ = E E (N x -x)(N y -\y\y ■ 1 



X=Nyy=-(Ny-l) 

x 2 — y 2 

y=l x =-(N y -l) 



( x 2 + y 2yo/2 



£ y (g2 + y2 y )5/2 (6) 



It is noteworthy that the previous expression is exact: its use for a numerical calculation 
of S requires the evaluation of 2N x N y terms, whilst Eq. (|J) demanded N 2 N 2 terms! 

The evaluation of S in the continuum approximation {Y^ x ,y ~* S J dxdy) is performed in 
the Appendix A. Here we will discuss the results. First of all, let us consider the case of 
an infinite aspect ratio (r = N x /N y = oo with N y finite). The quantity S per spin is (see 
Eq. (0)): 

An important feature of the previous expression immediately comes out: the shape 
anisotropy (per spin) vanishes upon increasing N y , even if N x /N y = oo! This means that 
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the shape anisotropy of an infinite stripe (N x = oo) becomes rapidly negligible, as N y 
increases. 

The numerical value of C\ cannot be determined by our "zero-order" continuum ap- 
proximation; in fact, the first-order correction, given by the term in square brackets of the 
Euler-Maclaurin summation formula (see Ref. 

£;/(*)= [ b dxf(x)+ 1 -[f(a) + f(b)} + ... (8) 

x=a Ja / 

contributes just to the term of order {1/N y ) in (|7J). Since the leading term (« (1/N y ) hxN y ) 
dominates only logarithmically, the constant cannot be neglected. In Eq. (|7|), c\ plays the 
role of the anisotropy for a single line of spins: 



S 

Cl 



N X Ny 



oo -i 

E73=C(3). (9) 
1=1 1 



So, we will rewrite (0) in the form: 



5 2 .niV s + #. (10) 



N x N y 3N y " y N % 

In the limit N x = oo, the numerical calculation of S/N x can be made much more 
efficient by means of the Ewald's summation technique [|IU|| , which allows to rewrite S/N x 
as an exponentially converging sum. This is done in Appendix B. We are therefore able to 
compare (see Fig. 2) the exact numerical result with the analytic approximation (|r0|) . Even 
our "zero-order" continuum approximation gives a fairly good approximation. 

The following considerations are meant to corroborate further on the previous results. 
If N x = oo, we can exploit the translational invariance in the x direction to write down: 

C 1 N v 00 /2 _ ( m _ m l\1 i N v 

— = - J2 [ ' = - V S(m - m!) 

N * 2 m,m'=l l=~oo [I 2 + (m- m') 2 ] 5/2 2 m ,m'=l 

= ~N y S(0) + (N y - c)S(c) 

Z c=l 



S(c) represents the "interaction" per unit length of two lines at distance c, whilst S(0)/2 

In the limit where N y 

d ( S \ S{0) 



is the self-interaction of a line. In the limit where N y can be treated as a continuous variable: 



c 



dN y \N X J 2 c=1 

If (I) ^ s ^ 

If the continuum approximation is applied also to the x direction, for c^Owe will have: 

/+oo £,2 2 

So, 

9 2 /5 



9iV2 y N j 3N 2 
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which gives 



— = ~\nN y + c N y + Cl 



(12) 



This expression would disagree with (0), if cq were not zero. It is zero, indeed. In fact: 



5(0) 



(13) 



c=l 



and by writing down 
system, we obtain: 



for a completely translationally invariant two dimensional 



N X Ny 



I 2 



m 



vJ^o) (^ 2 + - 2 ) 5/2 
S(0) 



m=l 



m) 



Co 



Therefore c vanishes and Eq. (pT2|) reduces to Eq. (|7|). 

After having discussed the shape anisotropy of an infinite stripe (N x = oo), upon increas- 
ing its "thickness" N y , now let us analyze the dependence of S/N x N y on the aspect ratio 
r. By handling the expression for „ S N given in Appendix and by using the "boundary" 



condition 



OO, Ny 



C(3), we obtain: 



N X Ny 



In N v 

3N, y 



03) 
N, 



F(r) 



(14) 



where F(oo) = 1 and F(l) « 1: the actual value of -F(l) is not really relevant, because the 
shape anisotropy vanishes in the limit r — 1. 

In Fig. 3 we compare the previous expression (as a function of r, keeping fixed the value 
of N y = 40) with the exact numerical result, obtained by exploiting Eq. (Q). The behavior 
is well reproduced by the analytical expression. 



3 Domain structures in a striped monolayer 

Now, let us turn to the study of domain structures in striped monolayers. Our purpose is 
to check if —and when— the appearance of magnetic domains is energetically favoured. We 
will consider a striped monolayer which is infinite in the x direction (N x = oo) and we ask 
for which values of the parameters, the creation of a domain wall along the y axis makes 
the system gain energy with respect to the collinear FM state. 

We will consider two cases, according to the direction of the magnetization: perpendicu- 
lar to the striped monolayer, and in the plane, along the "hard" direction (y). We consider 
just these two possibilities, because a striped monolayer has a "double" shape anisotropy: 
it has an easy-plane effect (as in a film, or in an infinite monolayer), but in the plane, it 
has also an easy- axis effect, along the x direction, which has been discussed in detail in the 
previous section. If the easy-plane effect is overcome by some anisotropics (K±) oriented 
along z, the film will break into domains. So, our first purpose (I) will be to study the finite 
size effects on this domain structure. 
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Conversely, if the easy-plane effect of Edi P prevails, the magnetization will be oriented 
along x, or —in presence of a sufficiently strong anisotropy K\\ favouring the y axis— along 
the smaller size of the stripe. In analogy with (I), we will expect that the striped monolayer 
breaks into in-plane domains. We will analyze (II) this possibility, by showing that domains 
don't appear, except for very large values of the dipolar coupling. 

What we have to do —in cases (I) and (II)— is to create a single domain wall along 
the y axis, and to compare the domain wall energy, with the dipolar energy gain. In the 
following, w will indicate the domain wall size. For an evaluation of the dipolar term, in the 
Heisenberg model it will be sufficient to consider Ising-like spins, with domains separated 
by an empty region of size w 



In both cases (out-of-plane and in-plane domains) the domain wall energy per unit length 
in the y direction, will be written as Edw (for the Ising model: Ed w = 2J, where J is the 
exchange coupling constant, and for the Heisenberg model: Edw — 2V JK, where K = K± 
in case (I) and K = K» in case (II)). 

I: out-of-plane domains. Because of the finite extension of the striped monolayer, the 
translational invariance in the y direction is lost and the dipolar interaction between a 
given spin and the spins of a neighbouring domain depends on the y coordinate of the 
spin. However, for the evaluation of the order of magnitude of the dipolar energy gain (per 
unit length), we can use the following approximate formula^, where — Q/[(x — x') 2 + y 2 ] 3 ^ 2 
is nothing but the dipolar interaction between two spins located in (x',0) and {x,y), and 
pointing along ±z: 



AE dip {N y ) 



-2Q 



-2Q 



In 



dx' 



dy 



dx 



m + w 2 



[{x - x'f + y 2 fl 2 



m + w 2 



w 



w 



N, 



(15) 



In the limit N y > w, AE, 

\AE dip \ > E du>} or 



dip 



-2ttln(2Ny/w). The FM 

'E dw ' 



w 



N y > — exp 



state will be destabilized if 



(16) 



2 ~" r V 2Q 

The explanation of this result is straightforward: out-of-plane domains will appear only 
if the lateral dimensions of the striped monolayer are larger than the size L that domains 
would have in an infinite monolayer. In fact, the right-hand side of Eq. (|16|) is exactly the 
typical size of a domain in an infinite monolayer Ul2 , l3[ . 

writes: 



II: in-plane domains, magnetized along ±y. In this case, AEdi P 

1 



AE d ip 



-2Q 



-2Q 



dx' 



dy / dx 



3y 2 



[(x 



x') 2 + y 2 } 3 / 2 



[(x 



x') 2 + y 2 fl 2 



w 



N, 



(17) 



where the dipolar interaction — Q |l/[(a; — x') 2 + y 2 ] 3 ^ 2 — 3y 2 /[(x — x') 2 + y 2 ] 5 ^ 2 | between 
two spins located in (x', 0) and (x, y) now has a contribution from the anisotropic part (the 
second term of Eq. ([[])). 

1 Indeed, it is possible to show that the dipolar energy gain is in between AE,np{Ny) and 2AE c i ip (N y /2). 
Since AE^ ip depends only logarithmically on 7V y , the correction is not important. 
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In the limit N y 3> w, Eq. simply writes: AE dip « — 2f2. So, the dipolar energy 
gain does not increase with N y and therefore the condition for the appearance of in-plane 
domains (\AE dip \ > E dw ) can be fulfilled only if Q « E dw , i.e. a rather large value. 

4 Conclusions 

Bulk systems differ substantially from two dimensional ones, because shape effects are 
relevant in three dimensions, but in-plane shape effects are negligible for stripes. More 
precisely, if the magnetization lies in the plane, its direction is mainly determined by existing 
magnetocrystalline anisotropies, because dipolar shape effects rapidly vanish when the sizes 
of the stripe increase, even if the aspect ratio goes to infinity. In particular, we have shown 
that the shape anisotropy per spin vanishes as ^j^- upon increasing the size N y , so that, 
e.g. for N y = 40 (and N x = oo) it is reduced by a factor larger than 10 with respect to a 
single chain of spins (N y = 1). 

The extreme weakness of shape effects in stripes has a further consequence on the ex- 
istence of in-plane domains: in fact, if the magnetization is forced in the y direction (the 
"hard" direction with respect to the shape anisotropy), in-plane domains with the magne- 
tization alternately directed along ±y appear only if Q is fairly large. 

We remark that our analysis has assumed that the striped magnetic monolayer consid- 
ered has a microscopic unit cell of square symmetry, so the theory is directly applicable 
to epitaxial monolayers grown on (100) substrates. In the case of the other high symme- 
try orientation - the (111) one - the overlayer has a triangular symmetry and no in-plane 
anisotropy is induced by the dipolar interaction, in the limit of an infinite monolayer. For 
striped monolayers, our previous treatment should be relevant. Conversely, for (110) sub- 
strates the microscopic unit cell has a rectangular symmetry, which induces a further (and 
possibly competitive) anisotropy in addition to the shape anisotropy. 

Acknowledgements - We acknowledge Danilo Pescia for having introduced us to the 
problem of the striped monolayers. 



A Evaluation of S in the continuum approximation 



In the continuum approximation, Eq. (0) for S rewrites: 



S 



2h - 2(N X - N y )I 2 , 




where: 



A(x,y) 



I 



i 




(x 2 + y 2 ) 5 / 2 
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The following integrals are easily calculated WM: 



A(x 1 ,x 2 ,yi,y 2 ) 



■1-2 



111 



dx I dyA(x, y) 

.rj Jy 1 

vl - A 



vl 



3y 1 x 2 Jyi + xi Sx 2 y 2 Jxi + yi 



[{xi,yi) ^ {x 2 ,y 2 )] (19) 



rX2 ry2 

B{x 1 ,x 2 ,yi,y 2 ) = dx dyxyA{x,y) 

Jx 1 Jy 1 



rp 2 



vl 



+ 



3^*2 + yi 



+ [(xi,yi) «-> (x 2 ,y 2 )} 



(20) 



C(x 1 ,x 2 ,y 1 ,y 2 ) 



X2 ry2 

dx / dy xA(x,y) 

xi Jyi 



3 



2/1 



2/i + ^2 



2/i 



2/i + z? 



+ 



22/1 



22/1 



+ [(^1,2/1)^(^2,2/2)] 



(21) 



D(xi,x 2 ,y u y 2 ) = dx dyyA(x,y) = -C(y 1 ,y 2 ,x 1 ,x 2 ) (22) 

► (^2,2/2)] means that we have to interchange x\ 



In the previous expressions, [(xi,yi) < 
with x 2 and y\ with y 2 - 

After some lengthy, but easy, calculations it is found that 



where 



and 



At 
B 1 
Ci 



h = N x N y A\ + B 1 - N x Di - N y d 
N 2 X - N 2 y V2 



3N X N V JN* + N* 3iV y 2 



N 2 x -N 2 y ^yg N y -N x 
3 3 



3 JiV2 + ^ 



1 



In 



I In 

3 



N y + JN* + Nl 



2N„ 



V2 

3 3 V / iVfTiV| 



iV, + ^ + N*) 
N x (l + v/2) 



V2 



+ 



2JV, 



3 3JW + N* 



h = - In 

3 



2iV, 



v/2 



As a function of the previous quantities, the shape anisotropy per spin writes: 



S 



h 



N„. — N, 



2N x N y 



N x N y 



N x N y 



y I 2 = A 1 + 



Bt D x d 



rNy 2 N y rN y 



1 - 



1\ h 



r) N„ 



(23) 



(24) 



(25) 
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In the limit N x — > oo it is immediately found that 



N X Ny 



(26) 



Conversely, for any value of r = N x /N y we obtain the following expression: 



S 



N X Ny 



1 - 



2 h,W. + ™ 



3iV, 



(27) 



where F'(r) is a function which depends only very weakly on the aspect ratio r: F'(l) ~ 
F'(oo) m 1. Indeed, the term proportional to F'(r) corresponds to the term of order 
(1/Ny) in (|26|): i.e. a term which cannot be determined consistently by the used "zero- 
order" continuum approximation. Imposing the "boundary condition" 

C(3), corresponds to put F'{r) = ((3)F(r), with F(oo) = 1 and F(l) ~ 1. 



B Exact calculation of S for N x = oo 



Let us consider the case of a stripe with an infinite aspect ratio r = N x /N y = oo, i.e. with 
N x — oo and finite iV^,. Since in this limit translation invariance is restored along the x 
direction, it turns out that the summation in Eq. (f|) can be rewritten 



S_ 

N~ x 



E E 



[m — m 



l\2 



2 m,m'=ii=-oo + (m - m') 2 ] 5/2 



-iV^(O) +J2(N y - c)S(c 



(28) 



c=l 



where S(c), the interaction per unit length of two lines at a distance c, is 



oo j2 ^2 oo 

^ (c) = £o a 5 + c2 ) 5/2 = tSo (p + c2 ) 3/2 



_ 2 V - 

zi^o (^ 2 + c 2 ) 5/2 



(29) 



The self-interaction of a line, S'(O), is readily evaluated in terms of the Riemann's zeta 
function ((x) 

OO oo 

S(0)= Yl T7i3 = 2 Ei = 2 C( 3 ) = 2- 1.202057 (30) 



/■ - x I' I" i=l ^ 



To evaluate S(c) for c > we use a method which was developed by Ewald |T(J for 
converting two-dimensional dipole sums to a rapidly converging form. First we take into 
account the identity [[14] : 

1 1 r°° 

/ ,. , (31) 



r (u) 



dt t v ~ x e~ at 



to rewrite S(c) as follows 



5(c) 



r(3/2) £z 
Next we employ the identity 



I /*oo I POO 

Wj^Jo r(5/2)^io 



-{l 2 +c 2 )t 



(32) 



E * 

i=— oo 



-IH 



'1/2 



E e 

n=— oo 



-(7rn) 2 /t 



(33) 
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so that 



S( c ) = Y dt e~ ch e' {nn)2/t 
Jo 



K 2c 2 ^_ t 



r(3/2) r(5/2) 



(34) 



The integrals can be exactly evaluated in terms of the modified Bessel functions K v (x) [14 

v/2 



r°° , fh\ v ' 

J dt e~ at e- b/t t v ~ x = 2 f - 1 tf„(2 



a6) (a, b > 0) (35) 



At last we obtain, for c > 1: 



Ait °° 1fi7r 2 00 

5(c) = — 2 M ^i(2vrc|n|)-^- ]T n 2 K 2 {27cc\n\) (36) 



n=— oo n.=— oo 



At this point we take into account the x — > expansion || K v (x) ~ |T(i/)(|x) y , so that 
the n = terms in the summation are found to give the finite contributions 

lim nKi(2ircn) 



n->o 2nc 
lim r? KiQ'K cn) 



n^o 2tt 2 c 2 
and finally we obtain for c > 1 

r) Q_ OO oo 

5(c) = -— + — £ n Xi(27rcn) - ^— £ n 2 K 2 (2ncn) (37) 

^ C C n=l ^ n=l 

The first term on the r.h.s. gives the main contribution, which coincides with the result of the 
continuum limit. The two n— summations converge very rapidly because the Bessel functions 
present an exponential decay for high values of their argument: K v (x) ~ \f^ e ~ x [^-~^'i^(^ h ' 2 ~ 
1) + O(^)] for x — > 00. In practice, excellent convergence is obtained summing about ten 
terms. 
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Figure 1: (a) Striped monolayer of size N x x N y . (b) Three-dimensional slab of size N x x 
N y x N z (with N z < N x ,N y ). 



Figure 2: Shape anisotropy per spin, S/ (N x N y ), for a striped monolayer with infinite aspect 
ratio (r = N x /N y = oo), as a function of iV r Dots: exact numerical results, derived from 



Eqs. (|28|) and ([37]). Line: analytical approximation, Eq. (10). 



Figure 3: Shape anisotropy per spin, S/(N x N y ), for a striped monolayer with fixed N y = 40, 
as a function of the aspect ratio r = N x /N y . Dots: exact numerical results, Eq. (^|). 
Dashed line: analytical approximation, Eq. (|T4T ) with F(r) = 1. Full line: asymptotic value 
(0.09153) of Eq. QTJ) for r -> oo. 
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